These are the versions of codes and data that reflect what is in the 4pi paper, Leka+
2025 (submitted).

The parameter-generating code from the AFT maps is a very large code bundle
that relies on SolarSoft IDL that produces numerous files.  Rather than 
provide the full infrastructure and all supporting files, we provide
IDL 'save' files that contain structures with all the parametrizations used in the paper,
and the population assignments (effectively the event lists).
The linking between AFT-ID (the unique identifier) and the AFT-map filenames
is also included in this save file.  In practice, this save file was created
once and then used repeatedly (as one can see from the codes) for the analysis.

In many instances in these codes, the "method" refers to an identifying sub-directory
and all other filenames are parallel to ease the analysis.  Thus, be mindful 
of overwriting if you do not preserve the directory structure at some level.

Content:

README                
CODE:
skill_score_guts.pro : computes skill scores from a probability array and matching event array
skill_score_plot.pro : plotting routine for Reliability plots
AFT_AR_to_FULL_DISK_ANALYSIS.pro  : driver for FD skill score analysis (except ROCs)
AFT_AR_ANALYSIS.pro  : driver for Region-based  skill score analysis (except ROCs)
make_roc_output.pro : computes arrays needed to make ROC plots, computes ROCSS
PAPER_plot_rocs.pro  : driver to make ROC plots (region)
PAPER_plot_FD_rocs.pro : driver to make ROC plots (FD)
PAPER_make_eg_plot.pro : Make Figure 1 

DATA:
AFT_HMI_GONG.tar
AFT_HMI_Baseline.tar              
Under each, parallel filenames and their contents:
ardata_tracked_2010_3x_NWRA.json : AFT AR-ID tracking file
[C,M]10min_Z00max_visMASTER_lat0.0hr_val24.0hr_nwra_aft_mag.sav: IDL save file with parameters used
in the analysis, and event list ("pop") for C+ and M+ respectively.
[C,M]10min_Z00max_vis[1,2,12]_lat0.0hr_val24.0hr_FLUX_TOT_NL_FLUX_TOT.nwra_aft_mag_npda: Files containing the 
probabilities from the NCI 2-Variable NPDE computed for the vis=1, vis=2, and vis=[1+2] AFT-AR's.
noaa_events_nwra_combined_curated_19960823_20240814.json: full cleaned / curated NOAA AR list
final_noaa_regions_nwra_combined_fixed_19960823_20240814.json: full cleaned  curated NOAA event list

